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Abstract 

Various series expansions have been developed for the two-layer, S = |, 
square lattice Heisenberg antiferromagnet. High temperature expansions are 
used to calculate the temperature dependence of the susceptibility and spe- 
cific heat. At T = 0, Ising expansions are used to study the properties of 
the Neel-ordered phase, while dimer expansions are used to calculate the 
ground-state properties and excitation spectra of the magnetically disordered 
phase. The antiferromagnetic order-disorder transition point is determined to 
be {J2/J\)c = 2.537(5). Quantities computed include the staggered magne- 
tization, the susceptibility, the triplet spin-wave excitation spectra, the spin- 
wave velocity, and the spin-wave stiffness. We also estimates that the ratio of 
the intra- and inter-layer exchange constants to be J2/J1 — 0.07 for cuprate 
superconductor YBa2Cu306.2- 

PACS Indices: 75.10.-b., 75.10J., 75.40.Gb 
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I. INTRODUCTION 



The magnetic properties of low dimensional systems have been the subject of intense 
theoretical and experimental research in recent years. It is by now well established that 
one-dimensional Heisenberg antiferromagnets with integer spin have a gap in the excitation 
spectrum, whereas those with half-integer spin have gapless excitations. The former have 
a finite correlation length, while for the latter it is infinite with the spin-spin correlation 
function decaying to zero as a power law. In 2-dimensions, the unfrustrated square-lattice 
Heisenberg model has long range Neel order in the ground state. It has gapless Goldstone 
modes as expected. In recent years much interest has focussed on systems with intermediate 
dimensionality and on questions of crossovers between d — 1 and d = 2, and the investi- 
gations showed that the crossover from the single Heisenberg chain to the two-dimensional 
antiferromagnet square lattice, obtained by assembling chains to form a "ladder" of increas- 
ing width, is far from smooth. Heisenberg ladders with an even number of legs (chain) show 
a completely different behavior than odd-leg-ladders. While even-leg-ladders have a spin 
gap and short range correlation, odd-leg-ladders have no gap and power-law correlations. 

This paper is concerned with a system between 2 and 3 dimensions: the S = ~ two-layer 
Heisenberg antiferromagnet, where each layer is composed of a nearest-neighbor Heisenberg 
model on a square lattice, and there is a further antiferromagnetic coupling between corre- 
sponding sites of each layer. The system can be described by the following Hamiltonian: 

H — J i y j ^ S a j ■ S a j + J2 ^ Si ; j • S 2j j (1) 

a=l,2 (ij) i 

where S a ,i is a S=l/2 spin operator at site % of the ath layer, and denotes a pair 

of nearest neighbor sites on a square lattice. Here J\ is the interaction between nearest 
neighbor spins in one plane, and J2 is the interaction between adjacent spins from different 
layers. We denote the ratio of couplings as y, that is, y = J 2 / ' J\- 111 the present paper, 
we study only the case of antiferromagnetic coupling, that is, both Ji and J 2 are positive, 
although some of the series are applicable to ferromagnetic coupling case. In the small 
J2/J1 limit, the model describes two weakly interacting 2D Heisenberg antiferromagnets, 
where each of them is Neel ordered and possesses gapless Goldstone excitations. In the large 
J 2 /Ji limit, pairs of adjacent spins from different layers form spin singlets separated from 
triplet states by a gap, ~ J 2 , and the system has no long-range order. Thus one expects 
there is a transition from ordered antiferromagnet to a quantum disordered spin liquid 
at a certain critical ratio of J2/J1. The bilayer Heisenberg antiferromagnets has recently 
attracted attention for several reasons. On the experimental side, the model is relevant 
to the understanding of the magnetic properties of cuprate superconductors YBa 2 Cu 3 6+x 
containing two weakly coupled Cu0 2 layers On the theoretical side, it is probably 

the simplest unfrustrated spin system which displays a quantum disordering transition of 
0(3) universality class, and Quantum Monte Carlo (QMC) simulation on this model do 
not suffer from a sign problem, making it possible to obtain precise data on the critical 
properties of the quantum transition. Hida |4| has performed a dimer expansion up to 6th 
order in Ji/J 2 , and he found the critical ratio to be (J 2 / ' Ji) c = 2.56, Gelfand || has recently 
extended the dimer expansion to 8th order, and obtains the critical ratio (J 2 /Ji) c = 2.54(2). 
Sandvik et al. [|]-|8| have carried out a detailed QMC simulations for the finite-temperature 
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properties of small (J2/J1) and near-critical (J2/ 1 Ji) c systems, and found a similar critical 
ratio: {J2/Ji) c — 2.51(2). On the other hand, a Schwinger boson mean-field calculation 
by Millis and Monien || resulted in a very large critical ratio: (J 2 / ' J\) c = 4.48; a self- 
consistent spin-wave theory [|l(J also predicts a large critical ratio: ( J 2 / J\) c = 4.3. Chubukov 
and Morr ||10| have argued that the discrepancies between the spin-wave results and the 
numerical simulations (or series expansions) have a physical origin and are related to the fact 
that in the spin-wave approach one neglects longitudinal spin fluctuations, and after taking 
longitudinal spin fluctuations into account, they got the critical ratio: ( J 2 / Ji) c — 2.73, which 
is reasonably close to the numerical results. The critical ratio obtained by the Schwinger- 
boson Gutzwiller-projection method JTTJ] is (J 2 / ' J\)c — 3.51. 

We have carried out extensive series studies of this system, via Ising expansions and dimer 
expansions at T = 0, and also by high temperature series expansions. Our results confirm 
the existence of a order- disorder transition at critical coupling ratio J2/J1 = 2.537(5). The 
complete spin-wave excitation spectra, spin-wave velocity, and some order parameters are 
computed. We also compare our results with previous calculations. 

This paper is organized as follows: In Sec. II the dimer expansions are presented, and the 
antiferromagnetic order- disorder transition points is determined precisely. In Sec. Ill, the 
Heisenberg-Ising model is studied via the Ising expansions, which allow for estimates of the 
domain stability of Neel-order phase, and also allow a direct comparison with experimental 
data of YBa 2 Cu 3 6+x . In Sec. IV we present the results of the high-temperature expansions, 
and our conclusion is presented in the last section. 



II. DIMER EXPANSIONS 

In the limit that the exchange coupling along the rungs J 2 is much larger than the 
coupling Ji within the plane, that is y 3> 1, the rungs interact only weakly with each other, 
and the dominant configuration in the ground state is the product state with the spin on 
each rung forming a spin singlet. The Hamiltonian in Eq. (JI]) can be rewritten as, 

H/J 2 = H + (l/y)V (2) 

where 

Hq = }^ Si t i ■ S 2 j 

i 

(3) 

We can construct expansion in 1/yby treating the operator Hq as the unperturbed Hamilto- 
nian and the operator V is treated as a perturbation. The linked-cluster expansion method 
has been previously reviewed in several articles [12|-[rjJ, and will not be repeated here. 

We have carried out the dimer expansions for the T = ground state energy per site 
Eq/N, the antiferromagnetic susceptibility Xi an d the lowest lying triplet excitation spectrum 
A(k x ,k y ) (odd parity under interchange of the planes) up to order The resulting 

power series in 1/y for the ground state energy per site Eq/N and the antiferromagnetic 
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susceptibility x are presented in Table I. Table of series coefficients for the triplet excitation 
spectrum A(k x ,k y ) would require inordinate amount of space to reproduce in print, it is 
available from the authors on request. Instead, we also present in Table I the series for the 
minimum energy gap m = A(jr,ir). The dimer expansion was firstly carried out by Hida 
||] up to order (1/y) 6 in 1992, and extended to order (1/y) 8 by Gelfand || recently Our 
results agree with these previous results, and extend the series by three terms. A list of 
19093 linked clusters of up to 12 sites contribute to the triplet excitation spectrum. 

To determine the critical point y c , we construct the Dlog Pade approximants |15 to the 



X and m series, and since the transition should lie in the universality class of the classical 
d = 3 Heisenberg model (our analysis also supports this), we expect that the critical index 
for x and m should be approximately 1.40 and 0.71, respectively. The exponent-biased Dlog 
Pade approximants || gives the critical point l/y c = 0.3942(8), or y c = 2.537(5). 

The spectra for some particular values of J\/ J2 are illustrated in Fig. 1, where the direct 
sum to the series at y = y c is indeed consistent with the Pade approximants that one can 
construct. To compute the critical spin-wave velocity v, one expands the spectrum A in the 
vicinity of wave-vector (71*, 7r) up to k 2 : 

A(7r _ k x , tt - k y )[l/y] = C(l/y) + D(l/y)(k 2 x + %) + ■■■ (4) 

and it is easy to prove || that the critical spin-wave velocity is equal to (2CD) 1 / 2 at y c . 
The series coefficients for 2CD in J1/J2 are listed in Table I. Using Pade approximants, 
one can estimate the properties at critical point y c : the spin-wave velocity is estimated to 
be [v I J 2 ) 2 = 0.559(2), or v/J 2 = 0.748(2), the ground-state energy per site is E /N = 
—0.4456(4) J 2 , the optical spin-wave gap (the gap at T = (0,0)) is A opt = 1.737(1) J 2 , and 
the gap at X = (ir, 0) is Ax = 1.368(8) J 2 , where the uncertainly is largely associated with 
the uncertainly in y c . 



III. ISING EXPANSIONS 

To construct an expansion about the Ising limit for this system, one has to introduce an 
anisotropy parameter x, and writes the Hamiltonian for Heisenberg-Ising model as: 

tf/Ji = H + xV (5) 

where 

a=l,2 (i J\ i a,i 

V = i^a,i^a,j + Sa,i$a,j) + V /X^l.i^2.i + ^l,i^2,i) ~ * X! € a,i^a,i (6) 

a=l,2 i a,i 

and e aj j = ±1 on the two sublattices. The last term in both H Q and V is a local staggered field 
term, which can be included to improve convergence. The limits x = and x — 1 correspond 
to the Ising model and the isotropic Heisenberg model respectively. The operator Ho is taken 
as the unperturbed Hamiltonian, with the unperturbed ground state being the usual Neel 
state. The operator V is treated as a perturbation. It flips a pair of spins on neighbouring 
sites. 
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The Ising series have been calculated for the ground state energy per site Eq/N, the stag- 
gered magnetization M, the uniform perpendicular susceptibility x±i an d the lowest lying 
triplet excitation spectrum A(k x , k y , k z ) for several ratio of couplings and (simultaneously) 
for several values of t up to order x 10 (the series for uniform perpendicular susceptibility x± 
is one order less). The resulting series for y = 0.5, 1, 1.5, 2.5 and t — are listed in Tables 
II and III (the series coefficients for triplet excitation spectrum are listed only for case of 
y = 1), the series for other value of y are available on request. 

We can also compute the spin- wave velocity v through a similar expansion of A{k x , k y , 0) 
in the vicinity of k x = k y = as in Eq.([|), and the series coefficients for 2CD in x are also 
listed in Table II. With this series and the series for the perpendicular susceptibility 
one can get the series for the spin-wave stiffness p s through the following "hydrodynamic" 
relation: 

Ps = v 2 x± ■ (7) 

There are two branches of the spin-wave dispersion. From the series one can see that 
the symmetric (optical) excitation spectrum A{k x , k y , 0) (that is, with even parity under 
interchange of the planes) is related to the antisymmetric (acoustic) excitation spectrum 
A{k x , k y , 7r) (that is, with odd parity under interchange of the planes) by 

A(k x , k y ,0) = A(tt - k y ,7i - k x ,7r) (8) 

so we only consider the antisymmetric excitation spectrum here. 

The conventional linear spin-wave theory |T|,|2]JT6[1 predicts the magnon dispersions as: 



A(k x , k y , tt) = 4SJ![(1 - il) + y(l + 7fc )/2] 1/2 (9) 

where 7^ = |[cos(fc x ) + cos(k y )]. Therefore one can see that the acoustic spin-wave gap 
A aco = A (71", 7r, 7i") is zero even for nonzero J2, the optical gap is 



A opt = A(0,0,tt) =ASp 1 J 2 (10) 

which is non-zero as long as J2 7^ 0, and the gap at X = (jr, 0, it) is 

A x = 4SJ 1 (l + y/2) 1 / 2 . (11) 

In a recent paper, Millis and Monien got the following relation for the optical spin- 
wave gap by using the Schwinger boson method: 

A 2 

j 2 = 2E! (12) 

16p s + 2A opt /7r 1 ; 



and so at the limit of J2 — > 0, one results that A opt = 4Sy JiJ 2 Z p , which differs from Eq.(|10D 

by a factor of J~Z^, where Z p = p s /(JiS 2 ) ~ 0.72. 

Having obtained the series we first attempt to identify critical points and to determine 
the nature of the phase diagram in the (x — y) plane. Naive Dlog Pade analysis reveals 
lines of singularities, as shown in Fig. 2. These are obtained consistently for the staggered 
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magnetization M and the uniform perpendicular susceptibility x- We note, in particular, 
that for the critical point x c is about equal to 1 at J2/J1 = 2.55, this is consistent with the 
critical point determined by the dimer expansions. 

At the next stage of the analysis, we try to extrapolate the series to the isotropic point 
(x = 1) for those values of the exchange coupling parameters which lie within the Neel- 
ordered phase at x — 1. For this purpose, we first transform the series to a new variable 

5 = l-(l-x) 1/2 (13) 

to remove the singularity at x — 1 predicted by the spin-wave theory. This was first proposed 
by Huse [|l7j] and was also used in our earlier work on the square lattice case [fTER . We then 



use both the integrated first-order inhomogeneous differential approximants [[OJ and the 
naive Pade approximants to extrapolate the series to the isotropic point 5 — 1 (x — 1). 
The results for the ground state energy per site Eq/N, the staggered magnetization M, 
the uniform perpendicular susceptibility x±> an d spin- wave velocity v , and the spin- wave 
stiffness p s are shown in Figs. 3-5. The results for J2 = presented in the figures are 
taken from our earlier work on the single-layer square lattice [fL8yi9| . The estimates of the 
ground state energy per site E /N from dimer expansions is also shown in Fig. 3. Note that 
in the classical spin- wave (large S) level, the uniform perpendicular susceptibility x±> the 
spin- wave velocity v, and spin- wave stiffness p s are found to be: 

x d—«a = _L (1 _ 1// 4) 

^classical = 2^25^(1 + V / A) 1 ' 2 (14) 
^classical = _ y 2^^ 

so in order to remove the trivial dependent of these quantities on y, Figs. 4-5 have been 
presented in terms of quantum renormalizations of the classical spin-wave values. We note 
that M, and p s first increase for small J2/J1, pass through a maximum at about y ~ 0.5, 
and then decrease for larger value of J2/J1, and vanishes at about J2/J1 = 2.55, which is 
consistent with the critical point determined by the dimer expansion. The reason that in the 
case of small J2/J1, the interlayer coupling enhances the antiferromagnetic long-range order 
is the system acquires a weak three dimensionality and quantum fluctuation is suppressed. 

The antisymmetric excitation spectra A(k x , k y , it) for some particular values of y are 
illustrated in Fig. 6, where we can see that the excitation is gapless at (tt, it, it) point. In 
Fig. 7, we present a plot of energy gaps at (0, 0, n) and (tt, 0, 7r) obtained from the linear 
spin- wave expansion, the Schwinger boson method 0, the dimer expansions, and the Ising 
expansions. We can see, from Figs. 6-7, that the spin-wave gap evolves very smoothly from 
the ordered side to the disordered side of the transition, (the gap at X = (0, 0, tt) turns 
out to have a small discontinuity at transition point, which may be due to the estimates 
from Ising expansions being a little too small). From Fig. 7, we can also see that for Ax, 
the results of leading order spin-wave theory is about 10% smaller than our estimates from 
Ising expansions. For optical magnon gap A opt , the Schwinger boson method give consistent 
results with those estimated from the Ising expansions over the region 0.2 ^ y ^ 1, but as y 
decreases from 0.2, the optical magnon gap estimated by the Ising expansions becomes larger 
than that predicted by the spin-wave theory, and also by the Schwinger boson method. 
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Finally, let us estimate the interlayer exchange constant for YBa2Cu306.2 using the opti- 
cal magnon gap estimated by Ising expansions: the recent neutron scattering measurements 
]T|fJ] have report that the optical magnon gap for YBa 2 Cu30 6 . 2 is A opt ~ 70meV, and if we 
take the canonical value of J\ = 120meV, we estimate that J2/J1 — 0.07, whereas J2/J1 
estimated by the spin-wave theory Jl],[2[ and by the Schwinger boson theory || are 0.08 and 
0.14, respectively. 



IV. HIGH TEMPERATURE SERIES EXPANSIONS 

We now turn to the thermodynamic properties of the system at finite temperatures. We 
have developed high-temperature series expansions for the uniform magnetic susceptibility 
x(T) and the specific heat C(T), for the system with arbitrary J\j J 2 , 



m 1 TrS?S^ H 

'\(n T H — 



(15) 



dU 
df 



C{T) 

where N is the number of sites, and (3 = J 2 /(kBT), and the internal energy U is defined by 



TrHe~ pH , s 

The series were computed to order (3 9 for x(T) and to order /3 8 for C{T). These series are 
listed in Table IV. We use the integrated first-order inhomogeneous differential approximants 



T5| and the naive Pade approximants to extrapolate the series. The resulting estimates for 
some particular values of Ji/J 2 ar e shown in Figs. 8-9. For the susceptibility, in order to 
remove the trivial dependence of x{T) on Curie- Weiss temperature, we present in Fig. 8 the 
(l/x(T) — T cw )/ Ji as a function of T j J\ where the Curie- Weiss temperature T cw = J x + J 2 /4. 
From this graph, we can see that for y > y c , the susceptibility x(T) converge nicely to 
very low temperature, and vanishes in the low temperature limit, which consistent with 
a gapped system. For y < y c , the extrapolation, as judged by the consistency of different 
approximants, do not converge well in low-temperature region. At the critical ratio {J2/ J\)ci 
Monte Carlo simulations have shown that the universal, linear behavior of x{T) extends 
up to very high temperature T ~ J±, unfortunately our series extrapolation is not reliable 
enough to see this behavior. 



V. CONCLUSIONS 



In the preceding sections we have presented and analyzed a variety of high-order pertur- 
bation series expansions for the two-layer, S — ^, square lattice Heisenberg (and Heisenberg- 
Ising) model. At T = 0, we have added three more terms to existing series for the ground 
state energy E /N, the antiferromagnetic susceptibility Xi an d the lowest lying triplet exci- 
tation spectrum A(k x ,k y ) from the expansions about the interlayer dimer limits, and also 
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carried out new perturbation series expansions about the Ising limit up to tenth order for 
a number of properties. In addition, high-temperature series for the system with arbitrary 
Ji/J 2 are computed to order /3 9 for uniform magnetic susceptibility x(T) and to order j3 8 
for the specific heat C{T). 

Based on exponent-biased analyses of the high-order dimer expansions for triplet energy 
gap and antiferromagnetic susceptibility, we estimate the antiferromagnetic order- disorder 
transition occurring at T = at a critical coupling ratio (J 2 / Ji) c — 2.537(5), which is 
in good agreement with the quantum Monte Carlo results || (J2/ Ji) c — 2.51(2). The 
spin-wave velocity at the critical point is estimated to be v/J 2 = 0.748(2), (equivalently, 
v/Jt = 1.898(9)), which is slightly higher than the QMC results § vjJ x = 1.69(2). The 
spin-triplet excitation spectra for the system on both ordered and disordered phases are 
computed, and they seem to be evolve smoothly from the ordered side to the disordered 
side of the transition. At the limit of Jil J\ — > 0, the optical magnon gap obtained by the 
Ising expansions is larger than 2^ JiJ 2 , instead of less than 2^ J\J 2 as predicted by the 
Schwinger boson method, and this implies the interlayer exchange constant is J 2 — 0.07Ji 
for YBa 2 Cu30 6 . 2 . 
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FIGURES 



FIG. 1. Plot of the antisymmetric spin-triplet excitation spectrum A(k x ,k y ) along 
high-symmetry cuts through the Brillouin zone for the system with coupling ratios 
J1/J2 = 0.1,0.2,0.3,0.3942 (shown in the figure from the top to the bottom at respec- 
tively), the lines are the estimates by direct sum to the dimer series, and the points (circles with 
error bar for the case of J1/J2 = 0.3942 only) are the estimates of the Pade approximants to the 
dimer series. 

FIG. 2. The phase boundary for the Heisenberg-Ising model without the staggered field (t = 0) 
obtained from analysis of the Ising series for staggered magnetization M and uniform perpendicular 
susceptibility x±- 

FIG. 3. The rescaled ground-state energy per site Eq/N as function of J2/ {Ji+Jz)- The crosses 
with dashed line connecting them are the estimates from Ising expansion, and the solid curve is the 
estimates derived from the dimer expansions. The error bars are much smaller than the symbols. 

FIG. 4. The staggered magnatization M and the renormalization factor of the uniform per- 
pendicular susceptibility Z x = 8xj_Ji/(1 — y/4) (at T = 0) versus J2/J1 as estimated by Ising 
expansions. 

FIG. 5. The renormalization factor of the spin-wave velocity Z c = v/[Ji(2 + y/2) 1 / 2 ] and the 
spin-wave stiffness Z p = 4p s /[Ji(l - y 2 /16)] versus J2/J1 as estimated by Ising expansions and 
dimer expansions (for Z c at the critical ratio (J2/ J\)c = 2.537 only). 

FIG. 6. Plot of the antisymmetric spin-triplet excitation spectrum A(k x ,k y ,ir) (derived from 
the Ising expansions) along high-symmetry cuts through the Brillouin zone for the system with 
coupling ratios J2/J1 = 0.5,1,1.5,2.5 (shown in the figure from the top to the bottom at (vr,0) 
respectively) 

FIG. 7. The rescaled energy gap at (0, 0, n) (the optical magnon gap) and at X = (ir, 0, n) as 
function of J2HJ1 + 'h)- The solid curves at large J2HJ1 + J2) are the extrapolation derived from 
the dimer expansions, the crosses with error bars are the estimates of Ising expansions, the long 
dashed lines at small J2HJ1 + J2) are the results of the linear spin-wave theory, the open circles 
connected by a short dashed line are the results of Eq. (|i~2^) (where the p s is taken from the results 
of Fig. 5), and the vertical dotted line indicates the position of the transition. 

FIG. 8. Susceptibility (l/x(T) — T cw )/J\ as a function of temperature for the system with 
coupling ratios J2/J1 = 5,10/3,2.537,5/3,1, several different integrated differential approximants 
to the high temperature series are shown for each J2/J1. 
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FIG. 9. The specific heat C(T) as a function of temperature for the system with coupling ratios 
J1/J2 = 0.1,0.3, several different integrated differential approximants to the high temperature 
series are shown for each J1/J2 
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TABLES 



■hx 



m/J2 



TABLE I. Series coefficients for dimer expansions of the ground-state ener; 
the antiferromagnet susceptibility x, the gap m/J2, and the critical spin- wave 
in the text). Coefficients of (l/y) n up to order n = 11 are listed. 

n 

~0~~ 

1 

2 
3 
4 
5 
6 
7 



E /(NJ 2 ) 
-3.750000000X10- 1 

0.000000000 
-3.750000000X10- 1 
-1.875000000X10- 1 
-1.171875000X10- 1 
1.289062500X10- 1 
2.373046875X10- 1 
2.162475586X10- 1 
-4.346008301 xlO" 1 
-9.492845535X10" 1 
-1.231647968 
6.411397507xl0 _1 



gy per site E /(NJ 2 ), 
velocity (as described 

2CD/J 2 
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10 
11 



1.000000000 
4.000000000 
1. 300000000 xlO 1 
3.950000000 xlO 1 
1.128125000X10 2 
3.131076389X10 2 
8.535167824 xlO 2 
2.304592044 x10 s 
6.169223926X10 3 
1.640000304 xlO 4 
4.332290008 xlO 4 
1.139206022x10 s 



1.000000000 
-2.000000000 
0.000000000 
-1.500000000 
-1.500000000 
-5.468750000X10- 1 
1.191406250X10- 1 
-4.202636719 
-1.060513306X10 1 
-2.693954849 xlO 1 
-3.315443579X10 1 
-8.044747781 xlO 1 



0.000000000 

1.000000000 

0.000000000 

2.375000000 

1.437500000 
-1.320312500 
-7.851562500X10- 1 
-2.716064453 

3.307617188 

1.409708977X10 1 

3.259252389 

2.037096158 xlO 1 
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TABLE II. Series coefficients for Ising expansions of the ground-state energy per site 
Eq/(NJi), the staggered magnetization M, uniform perpendicular susceptibility the energy 
gap m = A(0, 0, 0) = A(7r, it, it), and the critical spin-wave velocity 2CD (as described in the text) 
for J2/J1 = 0.5, 1, 1.5, 2.5 and t = 0, Nonzero coefficients x n up to order n = 10 are listed. 



11 


J2/J1 = 0.5 


Jo / 7, = 1 




J2/J1 = 1.5 


J0/J1 =25 







-5.62500000x10- 


1 


—6 25000000x10 


En/C 7i N) 
-1 


-6.87500000x10- 


1 


—8 12500000x10- 


1 


2 


-1.50669643x10- 


1 


— 1 56250000x10 


-1 


-1.81423611x10" 


1 


—2.86221591 x 10~ 


1 


4 


1.06782048x10- 


3 


1 3251 4881 x 1 0' 


-5 


-1.82797961x10- 


3 


— 1 01025331 x10~ 


2 


6 


-1.55432939x10- 


3 


-2.51031531x10 


-3 


-3.56180060x10- 


3 


-5.93087945x10- 


3 


8 


-6.14963212x10- 


4 


-8.46200313x10 


-4 


-1.01380621x10- 


3 


-1.82417072x10- 


3 


10 


-3.42105908x10- 


4 


-5.02769527x10 


-4 


-6.24621980x10- 


4 


-1.03444977x10- 


3 










JlX± 













2.22222222x10- 


1 


2.00000000x10 


-1 


1.81818182x10- 


1 


1.53846154x10- 


1 


1 


-2.81746032x10- 


1 


-2.50000000x10 


-1 


-2.29797980x10- 


1 


-2.08041958x10- 


1 


2 


2.92713075x10- 


1 


2.59090909x10 


-1 


2.37859757x10- 


1 


2.10928368x10- 


1 


3 


-3.04452465x10- 


1 


-2.68827110x10 


-1 


-2.49554660x10- 


1 


-2.34325217x10- 


1 


4 


3.05483508x10- 


1 


2.69863988x10 


-1 


2.49828011x10- 


1 


2.29571097x10- 


1 


5 


-3.09716015x10- 


1 


-2.74883091x10 


-1 


-2.56566684x10- 


1 


-2.43161906x10- 


1 


6 


3.10814223x10- 


1 


2.76075768x10 


-1 


2.57806919x10- 


1 


2.44310355x10- 


1 


7 


-3.13649358x10- 


1 


-2.79129287x10 


-1 


-2.61531350x10- 


1 


-2.53169737x10- 


1 


8 


3.14366380x10- 


1 


2.79764322x10 


-1 


2.62062184x10- 


1 


2.52501211x10- 


1 


9 


-3.16219779x10- 


1 


-2.81901242x10 


-1 


-2.64657763x10- 


1 


-2.58530073x10- 


1 






5.00000000x10- 


-1 


5.00000000x10 


M 

-1 


5.00000000 x 10" 


-1 


5.00000000x10- 


-1 


2 


-8.55389031x10" 


-2 


-7.81250000x10 


-2 


-8.45389660 xlO" 


-2 


-1.30714101x10- 


1 


4 


-8.46292613x10- 


-3 


-8.53938581x10 


-3 


-1.34925210x10" 


-2 


-4.35583219x10- 


2 


6 


-5.85842608x10- 


-3 


-8.07527219x10" 


-3 


-1.15477220x10- 


-2 


-2.85436734x10- 


2 


8 


-3.93110760x10- 


-3 


-4.86105744x10" 


-3 


-6.30805266 xlO" 


-3 


-1.93128872x10- 


2 


10 


-2.73750635x10- 


3 


-3.51123409x10 


-3 


-4.48806496x10" 


-3 


-1.37095158x10- 


-2 










m/Ji 













2.25000000 




2.50000000 




2.75000000 




3.25000000 




2 


-1.64330357 




-1.77083333 




-1.98338294 




-2.55445076 




4 


1.20240301x10- 


-1 


5.95393105x10 


-2 


9.76911256x10- 


-2 


3.91257865x10- 


1 


6 


-2.91150072x10- 


-1 


-2.85968590x10" 


-1 


-3.69547989x10- 


-1 


-9.18026495x10- 


■ 1 


8 


8.98232628x10- 


-2 


6.28691338x10 


-2 


1.56753567x10- 


-1 


1.10864903 




10 


-1.90221034x10" 


-1 


-1.73683175x10 


-1 


-3.02886458x10- 


-1 


-2.00750062 




2 


3.97500000 




4.16666667 


2CD/f( 


4.51785714 




5.59722222 




4 


-1.62276149 




-1.51750579 




-1.78959346 




-3.20894532 




6 
8 


1.30765116 
-1.40693583 




1.08647250 
-1.13102684 




1.39046692 
-1.70310223 




3.71475815 
-7.26139043 




10 


1.60579536 




1.23771527 




2.19444652 




1.47488833X10 1 
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TABLE III. Series coefficients for the Ising expansions of the triplet spin-wave excitation 
spectrum A(k x , k y , k z ) = J x Y,k,n,m, P a k,n,m, P x k cos(pk z )[cos(mk x ) cos(nk y ) + cos(nk x ) cos(mk y )]/2. 



Nonzero coefficients ak^ n ,m,p U P to order k = 10 for the case of J2 = J\ are listed. 



(k,n,m,p) 


^k.n^m^p 




(k,n,m,p) 






(k,n,m,p) 


^k ,n.m,p 




(k,n,m,p) 


,n.m,p 




( 0,0,0,0) 


2.50000000 




( 6,1,1,0) 


-5.01421733x10" 


■2 


( 6,3,3,0) 


-3.50360695x10" 


-3 


(10,7,0,1) 


-1.95483982x10" 


4 


( 2,0,0,0) 


-1.04166667x10" 


1 


( 8,1,1,0) 


3.13641370x10" 


-2 


( 8,3,3,0) 


-2.76765928x10" 


-3 


(10,9,0,1) 


-3.18431828x10" 


-6 


( 4,0,0,0) 


8.82430142x10" 


-2 


(10,1,1,0) 


-2.49468947x10" 


-2 


(10,3,3,0) 


-2.34255921x10" 


3 


( 4,2,1,1) 


-1.18055556x10" 


1 


( 6,0,0,0) 


-2.63720521x10" 


-2 


( 4,3,1,0) 


-3.93518519x10" 


-2 


( 8,5,3,0) 


-7.65220675x10" 


4 


( 6,2,1,1) 


-2.64155861x10" 


-2 


( 8,0,0,0) 


1.81429057x10" 


-2 


( 6,3,1,0) 


-1.87035451x10" 


-2 


(10,5,3,0) 


-1.48448940x10" 


3 


( 8,2,1,1) 


5.61540297x10" 


-3 


(10,0,0,0) 


-9.88646062x10" 


-3 


( 8,3,1,0) 


-3.19610735x10" 


-3 


(10,7,3,0) 


-7.64236388x10" 


5 


(10,2,1,1) 


-2.48059104x10" 


-2 


( 2,2,0,0) 


-3.33333333x10" 


1 


(10,3,1,0) 


-1.18183470x10" 


-2 


( 8,4,4,0) 


-4.78262922x10" 


-4 


( 6,4,1,1) 


-1.05108208x10" 


-2 


( 4,2,0,0) 


1.52199074x10" 


-2 


( 6,5,1,0) 


-2.10216417x10" 


-3 


(10,4,4,0) 


-8.91369509x10" 


4 


( 8,4,1,1) 


-5.39907287x10" 


-3 


( 6,2,0,0) 


-3.34161849x10" 


-2 


( 8,5,1,0) 


-2.49774309x10" 


-3 


(10,6,4,0) 


-1.33741368x10" 


4 


(10,4,1,1) 


-5.30710289x10" 


-3 


( 8,2,0,0) 


1.25625535x10" 


-2 


(10,5,1,0) 


-2.84123845x10" 


-3 


(10,5,5,0) 


-8.02448208x10" 


5 


( 8,6,1,1) 


-7.65220675x10" 


4 


(10,2,0,0) 


-1.75733942x10" 


-2 


( 8,7,1,0) 


-1.09317239x10" 


4 


( 2,1,0,1) 


-6.66666667x10" 


-1 


(10,6,1,1) 


-1.16719696x10" 


-3 


( 4,4,0,0) 


-4.91898148x10" 


.3 


(10,7,1,0) 


-3.13087047x10" 




( 4,1,0,1) 


6.74768519x10" 




(10,8,1,1) 


-5.73177291x10" 




( 6,4,0,0) 


-4.91156561x10" 


-3 


(10,9,1,0) 


-6.36863657x10" 


-6 


( 6,1,0,1) 


-5.84006701x10" 


-2 


( 6,3,2,1) 


-2.10216417x10" 


-2 


( 8,4,0,0) 


-2.30378946x10" 


-3 


( 4,2,2,0) 


-2.95138889x10" 


-2 


( 8,1,0,1) 


3.34318983x10" 


-2 


( 8,3,2,1) 


-5.04717657x10" 


-3 


(10,4,0,0) 


-3.17090252x10" 


-3 


( 6,2,2,0) 


-1.09825781x10" 


-2 


(10,1,0,1) 


-3.12876108x10" 


2 


(10,3,2,1) 


-7.50314475x10" 


-3 


( 6,6,0,0) 


-1.75180347x10" 


4 


( 8,2,2,0) 


-7.99539319x10" 


-4 


( 4,3,0,1) 


-1.96759259x10" 


-2 


( 8,5,2,1) 


-2.29566203x10" 


-3 


( 8,6,0,0) 


-3.97833036x10" 


4 


(10,2,2,0) 


-7.65794740x10" 


-3 


( 6,3,0,1) 


-1.30043279x10" 


-2 


(10,5,2,1) 


-2.08196519x10" 


-3 


(10,6,0,0) 


-6.37350643x10" 


4 


( 6,4,2,0) 


-5.25541042x10" 


-3 


( 8,3,0,1) 


-9.25242372x10" 


-4 


(10,7,2,1) 


-2.29270916x10" 


4 


( 8,8,0,0) 


-6.83232746x10" 


-6 


( 8,4,2,0) 


-4.57905354x10" 


-3 


(10,3,0,1) 


-6.65895455x10" 


3 


( 8,4,3,1) 


-3.82610338x10" 


-3 


(10,8,0,0) 


-3.37994352x10" 


-5 


(10,4,2,0) 


-4.19352456x10" 


-3 


( 6,5,0,1) 


-1.05108208x10" 


-3 


(10,4,3,1) 


-2.63082695x10" 


-3 


(10,10,0,0) 


-3.18431828x10" 


-7 


( 8,6,2,0) 


-3.82610338x10" 


-4 


( 8,5,0,1) 


-1.65065868x10" 


-3 


(10,6,3,1) 


-5.34965472x10" 


4 


( 2,1,1,0) 


-6.66666667x10" 


1 


(10,6,2,0) 


-8.44976702x10" 


-4 


(10,5,0,1) 


-1.45569445x10" 


3 


(10,5,4,1) 


-8.02448208x10" 


4 


( 4,1,1,0) 


1.00115741x10" 


1 


(10,8,2,0) 


-2.86588646x10" 


-5 


( 8,7,0,1) 


-5.46586197x10" 


-5 









TABLE IV. Series coefficients for the high-temperature expansions of the uniform suscepti- 
bility X (T) = T- 1 J2 min a min {J 1 /J 2 ) m P n and the specific heat C(T) = /3 2 £™,n a m , n (Ji/ J 2 ) m /3". 
Nonzero coefficients a m ,n up to order n = 9 for x or order n = 8 for C are listed. 



(m 


n) 


Q>m,n 




(m,n) 


Q>m.n 




(m, 


n) 






(m,n) 






( 0, 


0) 


2.500000000x10" 


1 


( 4, 4) 


1.692708333x10" 


X 
-2 


(T) 

( 0, 


7) 


-3.986661396x10" 


-5 


( 6, 8) 


1.061827644x10 


-2 


( °. 


1) 


-6.250000000x10" 


2 


( 0: 5) 


2.115885417x10" 


-4 


( h 


7) 


5.113389757x10" 


4 


( 7, 8) 


4.403056796x10 


-3 


( 1, 


1) 


-2.500000000x10" 


1 


( 1, 5) 


-3.580729167x10" 


-3 


( 2, 


7) 


-2.441406250x10" 


-3 


( 8, 8) 


3.978426494x10 


-4 


( °. 


2) 


-1.562500000x10" 


2 


( 2, 5) 


1.822916667x10" 


-2 


( 3, 


7) 


4.155815972x10" 


-3 


( 0, 9) 


3.899804709x10 


-6 


( 1, 


2) 


1.250000000x10" 


1 


( 3, 5) 


-2.864583333x10" 


-2 


( 4, 


7) 


6.287977431x10" 


-3 


(1,9) 


-4.881601485x10 


-5 


( 2, 


2) 


1.250000000x10" 


1 


( 4, 5) 


-4.492187500x10" 


-2 


( 5, 


7) 


-1.420898438x10" 


-2 


( 2, 9) 


2.027723524x10 


-4 


( 0, 


3) 


1.302083333x10" 


3 


( 5, 5) 


-9.244791667x10" 


-3 


( 6, 


7) 


-1.158582899x10" 


-2 


( 3, 9) 


-1.131895358x10 


-4 


( 1, 


3) 


1.562500000x10" 


2 


(0, 6) 


-1.044379340x10" 


-4 


( 7, 


7) 


6.285652282x10" 


4 


(4, 9) 


-1.836698017x10 


-3 


( 2, 


3) 


-1.250000000x10" 


1 


( 1, 6) 


5.533854167x10" 


-4 


( 0, 


8) 


1.513768756x10" 


-6 


( 5, 9) 


3.601243761x10 


-3 


( 3, 


3) 


-4.166666667x10" 


2 


( 2, 6) 


1.399739583x10" 


-3 


( 1 , 


8) 


3.700861855x10" 


-5 


( 6, 9) 


3.023339327x10 


-4 


( 0, 


4) 


1.627604167x10" 


3 


( 3, 6) 


-1.527777778x10" 


-2 


( 2 - 


8) 


-6.515260727x10" 


4 


( 7, 9) 


-7.574656653x10 


-3 


( 1, 


4) 


-1.041666667x10" 


2 


( 4, 6) 


2.438151042x10" 


-2 


( 3, 


8) 


3.424169147x10" 


-3 


( 8, 9) 


-1.439339774x10 


-3 


( 2, 


4) 


7.812500000x10" 


3 


( 5, 6) 


2.350260417x10" 


-2 


( 4, 


8) 


-5.407569522x10" 


-3 


( 9, 9) 


-5.211347415x10 


-4 


( 3, 


4) 


8.333333333x10" 


2 


( 6, 6) 


1.991102431x10" 


-3 


( 5, 


8) 


-1.251705109x10" 


-3 


















C(T) 














( 0, 


0) 


9.375000000x10" 


2 


( 2, 4) 


-1.025390625x10" 


-2 


( 3, 


6) 


3.450520833x10" 


-3 


( 9, 7) 


-6.530006045x10 


-2 


( 2, 


0) 


3.750000000x10" 


1 


( 3, 4) 


-1.367187500x10" 


-2 


( 4, 


6) 


-4.157307943x10" 


-2 


( 0, 8) 


-5.264736357x10 


-5 


( 0, 


1) 


4.687500000x10" 


2 


( 4, 4) 


1.464843750x10" 


-1 


( 5, 


6) 


2.871093750x10" 


-2 


( 2, 8) 


-5.870092483x10 


-4 


( 3, 


1) 


1.875000000x10" 


1 


( 6, 4) 


3.173828125x10" 


-2 


( 6, 


6) 


-8.115234375x10" 


-2 


( 3, 8) 


-4.486810593x10 


-4 


( 0, 


2) 


-5.859375000x10" 


3 


( 0: 5) 


9.399414063x10" 


-4 


( 8, 


6) 


1.951700846x10" 


-2 


( 4, 8) 


3.243491763x10 


-3 


( 2, 


2) 


-4.687500000x10" 


2 


( 2, 5) 


4.785156250x10" 


-3 


( 0, 


7) 


-1.816522507x10" 


-5 


( 5, 8) 


-1.445297968x10 


-2 


( 4, 


2) 


-1.640625000x10" 


1 


( 3, 5) 


3.759765625x10" 


-3 


( 2, 


7) 


7.376534598x10" 


4 


( 6, 8) 


2.286320641x10 


-2 


( 0, 


3) 


-9.765625000x10" 


3 


( 4, 5) 


4.990234375x10" 


-2 


( 3, 


7) 


2.947126116x10" 


4 


( 7, 8) 


-3.248014904x10 


-2 


( 2, 


3) 


-3.906250000x10" 


2 


( 5, 5) 


9.365234375x10" 


-2 


( 4, 


7) 


-2.781633650x10" 


-2 


( 8, 8) 


3.321502322x10 


-2 


( 3, 


3) 


-3.906250000x10" 


2 


( 7, 5) 


1.268066406x10" 


-1 


( 5, 


7) 


-2.221854074x10" 


-2 


(10, 8) 


-2.914089021x10 


-2 


( 5, 


3) 


-1.953125000x10" 


1 


( 0, 6) 


4.185994466x10" 


-4 


( 6, 


7) 


-2.977120536x10" 


-2 








( 0, 


4) 


-1.586914062x10" 


3 


( 2, 6) 


4.451497396x10" 


-3 


( 7, 


7) 


-1.089634487x10" 


1 
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